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I. INTRODUCTION 


Magnetic oscillations or de Haas-van Alphen effect (dHvA) in quasi-two-dimensional met¬ 
als are well accounted for by the Lifshitz-Kosevich (LK) theory [l-5|, which relates the 
frequency to the surface area enclosed by the cyclotronic trajectories. This geometrical 
approach is based on the semi-classical quantization theory of Onsager j^, and allows for 
the determination of many physical parameters of the Fermi surface (FS). Characteristic 
field Bf. = <Fo/Ao associated with the quantum flux trough the unit cell area Aq is generally 
very large. Indeed, for organic conductors with unit cell area as large as Aq = lOOA^, 
is still 4136 T. Therefore, available magnetic fields stay well within the limit of the semi- 
classical approximation. In the opposite case, e.g. for significantly larger unit cell or applied 
magnetic helds, quantum corrections to the Landau spectrum or modihcation of the 
Lifshitz-Kosevich (LK) theory would be necessary. 

A question that remains to be addressed deals with the effect of departure from parabolic 
curvature of the electronic band dispersion on the amplitude and phase of quantum oscil¬ 
lations. This question is particularly relevant in the case of Dirac fermions, the electronic 
dispersion of which is linear. In this case, the LK calculation based on a truncation at first 
order in energy when evaluating the grand potential is no more sufficient, since the Landau 
level energy is known to display a square root dependence on the Landau level index Q]. 

In this paper, we address the question of the relevance of the non-parabolicity in two 
classes of materials. First, the LK calculation for magnetization in the presence of a uniform 
field is reconsidered in the case of Landau quasiparticles, or conventional fermions, relevant 
to e.g. organic conductors. Then Dirac fermions, which have linear band dispersion, are 
considered. In both cases the field- and temperature-dependent phase factor of the Fourier 
coefficients is evaluated as a function of the FS curvature which is the main parameter of 
the model. We can also mention other work dealing with non-parabolicity of the Fermi 
surface, in the special case of a tight-binding model in two dimensions ^ where the band 
gap closes for certain hlling factor, and where the temperature amplitude does not follow 
the LK formula. 


The main results of this paper can be summarized as follow. For a band dispersion in 
a two dimensional material with a closed Fermi surface, we assume that the area of such 
surface is given by the quantity S{E) at energy close to the Fermi energy ~ p. The 
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effect of band curvature (defined as the second derivative S"{E) = d‘^S{E)/dE'^ near the 
Fermi energy E = fi) is mainly to add a phase shift contribution in each harmonics of 
order p of the oscillating quantities such as the magnetization. More precisely, the phase shift 
is a function of held and temperature and is given by the following semi-classical expression 
in the case of organic conductors 


— 

p{x) = 


Ti^ksTh^S {p) 

-(p 

m* 

sinh( 2 a;) — x — 



xcosh^(a:) /sinh^(x) 


( 1 ) 


where m* is the effective mass of the quasi-particle and Uc = eB/m* the cyclotronic 
frequency. In the case of Dirac fermions, for which the Landau energy level spectrum En 
increases like the square root of the index n, the expression of the phase shift is given instead 
by 


0p — 

ip{x) = 


71 ksT 
2 m*Vp 


(f [277 p- 


knT 


hUr 


2 sinh( 2 x) — x — xcosm(x) —2x^ sinh"^(x) / sinh"^(x) 


( 2 ) 


where vp is the Fermi velocity. In the following we will derive both Eq. ([T]) and Eq. ([2]) 
using semi-classical analysis, and study their asymptotic limits, when both T or B are varied 
within physical ranges. This effect has to be differentiated from other possible contributions 
coming from additional physical mechanisms. For example, the presence of a spin-orbit 
coupling leads to a splitting of the dHvA frequency whose magnitude is proportional to B‘^ 
and effective mass m* [^. Freq uency splitting due to spin-orbit coupling has been studied 
in bilayer high-T^ cuprates [ll| and can also be attributed to the splitting of the Fermi 
surface. Additional phases exist in the presence of magnetic breakdown. It is well known 
that magnetic breakdown is accompanied with a held-dependent phase |^, especially for 
large orbits. This phase depends more precisely on the ratio between the held and the 
breakdown held, and not temperature. Existence of this Onsager phase has been studied 
and revealed in organic compound 6 '-(BEDT-TTF) 4 CoBr 4 (C 6 H 4 Cl 2 ) where the breakdown 
held is close to 35T jl^. More recently, the existence of phase ohsets 7 has been questioned 


in t 


lermodynamic quantities of three-dimensional topological insulators with surface states 


The topological nature of these insulators can be detected due to the presence of 


a Berry phase within the oscillations when the particle-hole symmetry is broken and the 
material has a band gap [l^. These global phases (as for the Maslov index 7 , see below) 
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are independent of temperature, but present a linear variation with the held Q , contrary 
to the held and temperature-dependent phases Eq. ([1]) and Eq. ([2]) which come from a local 
ehect of the band dispersion as we will see in the next section. Non-zero topological Berry’s 
phase was also investigated in graphene 18| by measuring the magnetoresistance in the 
quantum Hall regime, with a compelling evidence of a value for 7 diherent from 1/2 due to 
the presence of Dirac fermions. Eor a review of the topological phases in condensed matter 
physics, see |^. 


II. EFFECT OF BAND CURVATURE CORRECTIONS IN ORGANIC 

CONDUCTORS 


Within the semi-classical framework, the phase quantization can be expressed in terms 
of surface swept by the quasiparticle in the Brillouin zone. It is given by the integral 

S{En = E) = ^(f kydk^ = ±b{n +j), (3) 

4vr^ JE=E(k^,ky) 

where b = eB/h is an ehective Planck constant or reduced held, and E^ is the energy of 
the Landau band with n the Landau level index. 7 is the Maslov index which is equal to 
1/2 for Landau quasiparticles with a parabolic band. In such a case, S{E) = m*E/{2TTh^) 
varies linearly with the energy. For massless Dirac electrons with a linear energy dispersion, 
S{E) = ttE‘^/( 2TThvF)‘^ and 7 = 0 . In two-dimensional systems, the grand potential is 
expressed by 

, 00 

n{B, p, T) = -A- log(l + exp[/3(/i - En)]), 

P n=0 

where /i is the chemical potential and A is the sample area. bA is the degeneracy of each 
Landau Level. The Poisson formula can be used for any discrete series over positive integers 
n: = /o°° E(? 7 ,) n -I- 2 /o°° F{n) cos(27rpn) n. This allows us to rewrite the 

oscillating part of the grand potential in terms of Fourier modes 

log(l exp[/3(p - En)]) exp{2i-Kpn) n. (4) 



Using exp( 2 i 7 rpn) —)• exp{2i7ipn)/{2i7ip) as primitive function, a double integration by parts 
can be performed, yielding 




osc 



4cosh^[|(E„ 


exp( 2 ^ 7 rpn) 

/i)] (2z7rp)2 


A 
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PKEi, 


ft 

n-^n 


exp(2i7rpn) 


4cosh^[f(E„-/i)] (2i7rp)= 


n. 


(5) 


where E'^ = dEn/dn and E'^ = d'^En/dv?. Since for a parabolic band i?" = 0, only the 
hrst term does not vanish in this case and E'^ = hooc is n-independent and proportional to 
b. In general E”, E”', etc... are non zero, and Eq. ([5]) can be solved aronnd E = fi nsing the 


formal series expansion (see also [^) 


ac 02 C 

S{E) = 5(/r) + {E- — + \{E - + • ■ • 


( 6 ) 


In the standard LK theory, only the first two terms in Eq. ([6]) are taken into acconnt: the 
hrst one typically sets the freqnency of the oscillations, while the second one (linear in the 
energy difference) gives rise to the thermal rednction factor (see page 184 of reference ^ 
after eqnation 10.28 for instance). The third term stndied in this manuscript, and quadratic 
in energy difference, is responsible for a phase shift of the oscillations as we will see further 
below. To illustrate the discussion, we can take a typical example, the tight binding model 
of free electrons on a discrete lattice with hopping parameter t and which is described by the 
energy dispersion E^k^, ky) = —t[cos{kx) + cos(/Cj^)]. The density of states can be written as 


AoS{E) = - 


TT^ 


-E/t 


'-2 


-K 


u 


2 + u 


u 


u, 


(7) 


where Aq is the area of the unit cell and K{k) is the complete elliptic integral of the hrst kind: 
K{k) = Jg u[(l — v?){l — Eq. ([7]) is computed from Eq. (E]) using the coordinate 

equation ky = cos“^(—E/t — cos(A;a;)), then differentiating it with respect to the energy to 
get S\E). After the change of variable u = cos(/cx), Eq. ([7]) can be rewritten as 


A, 


dS{E) 


u 


( 8 ) 


dE ttH ^(l-u^)[l-(u + E/ty] 

After the additional change of variable v = u + \E/t, the integral is symmetric around the 
origin, and a further transformation w = v/{l + \E/t) leads to an expression involving the 
elliptic integral. We can in particular perform an expansion around the lower band limit 
E = —2t such that 


^oS{E) —+ 2] + + 2) + 


1 /E 


(9) 


IbTT V t 

We can identify S'{E) = (27rAot)“^ with 2nm*/fi?, since At = Ah?/{AQim*) is the bandwidth. 
The curvature parameter of the surface area enclosed by the orbit dehned in the following 
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hj K = S"{E), can be rewritten as k = (SvrAot^)”^. Therefore for a given dispersion, we 
can relate the different coefficients of the expansion in Eq. ([ 6 ]) with microscopic parameters 
such as the hopping constants, effective mass or bandwidths. 

In the most general case, we would like to use the expansion Eq. (| 6 ]) up to second order 
to compute the different thermodynamic quantities such as Eq. (jll). Up to now, the LK 
calculation considers only the hrst order around the Fermi surface, whose coefficient is given 
by the effective mass (the slope of S{E)). The second order term will modify, as we will 
see below, mostly the phase of the magnetic oscillations. To include the effect of the second 
order term, we perform hrst a change of variable n —)• E in the expression of the grand 
potential Eq. ([5]) using dE = E'^dn 




A 


-2b Re 

1 


p=i 


/3K 


exp[ 2 z 7 rpn(E)] 


4cosh^[|(E -/r)] ( 2 z 7 rp)^ 


exp[ 2 z 7 rpn(E)] 


4cosh^[|(E -/r)] ( 2 z 7 rp)^ 


E. 


We then use the relations 


EL = 


tuVr 


E'' = 


n^tvjJcY 


2TTh^ 


K 


K = 


( 10 ) 


1 + k(E„-p)’ [1 + k(E„ - p)]2’ m" 

to obtain the oscillating part of the grand potential expressed as an integral over the energy 




A 


OO « 

-26ReV / 




exp \2inpn{E)] 


4[1 + k{E — /i)] cosh^[|(E — /i)] 

k/3{hjjJcY exp[ 2 i 7 rpn(E)] 

4[1 + k{E - /i)]2 cosh^[|(E - p)] (2i7rp)3 


{2i7Tpy 


E, 


( 11 ) 


where n{E) = S{E)/b — 7 , S{E) being given by the expansion Eq. ([ 6 ]) around the Fermi 
energy. The next step is to perform the integration around the saddle point E = fi at low 
temperature, using the variable x = (3{E — p)/2tt and replacing n{E) by n{x) 


, , S(u) 2 'kx 2'k‘^x‘^k 

n{x) = -7 + + 


( 12 ) 


b ' /3huc ' fd'^hojc 

In the case of a parabolic band, the hrst and second terms are related to the Onsager phase 
or Maslov index, and the oscillation frequency E = hS{fi)/e, respectively. More specihcally, 
magnetization is obtained by diherentiation of Eq. (ITTll with respect to minus B, yielding 


e^E ^ W 


Tim 


* / ^ 


^27rp 


sm 


2:rp— - 27r7p + 


(13) 
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where Ap and cpp are respectively the amplitude and phase of the imaginary damping factor 
Zp dehned by 


Zp = Apexp{i(f)p) 



1 -I- Tx/a + oox^ ja 

2mpx{T -f oox') 

2 J-j3fi/{2n) 

1 -|- ^ooxjr 

cosh^(7rx) 


and which involves the dimensionless parameters 


(14) 


F 27r TTTK 

B' f3hxJc /d 


(15) 


Hence, the temperature dependence is given both by r oc T and u (x while u, which 
is proportional to S”{E), includes the band curvature contribution. The integral can be 
computed by extending the lower bound to —oo at low temperature and by neglecting the 
factors in front of the exponential that are proportional to parameter which is generally 
small compared to unity for magnetic helds strength currently available, except for very small 
orbits area. Besides, the denominator (1 + 27rxK//3) is strictly positive as long as kju < 1. In 
the peculiar case where k = l/p, which is relevant for Dirac fermions, as developed below, 
we can show that u = r'^/Aa and that the overall factor reduces actually to 1 + Tx/{2a) 
with no divergence around the lower bound of integration. Following the LK theory and 
using the residue theorem, we can consider an integration path that goes over the upper 
complex plane when x is complex. Indeed, the singularities of the cosh(7ra;)“^ function give 
the main contributions to the integral since they are located on the positive imaginary plane 
Xn = i{n + ^). The oscillation amplitude can be written as an infinite summation over the 

Xn’s 


Zp = < 27rpr + Ainpu ( ^ j \ exp — 27rpr ( ^ j — 2i7ipu in + ^ 


(16) 


n=0 


For K = 0 (or a; = 0), the summation can be performed and we obtain the well known LK 
thermal reduction factor Zp = Rp = pirr/ sinh(p7rr) 20|, to which the Dingle factor and, 
eventually, magnetic breakdown probabilities are added. Eq. fITB]) can be expanded up to 
the hrst order in a; oc k, the sum over n performed, and the result re-exponentiated, which 
allows us to rewrite Eq. flTHl) as 


PTTT 

Zj — - 

^ sinh(p7rr) 


exp 


'. 00 

l — p['KpT 


= i?pexp 


i—nksTpi^To'^p- 


fvojr 


’p{x) = sinh(2a:) — x — xcosh^(a:) / sinh^(a;). 


( 17 ) 


















We can identify the amplitude Ap with Rp at this order. As a consequence, while the oscil¬ 
lation amplitude, which is accounted for by the real part of Zp, is unaffected by deviations 
from parabolicity, an extra phase is present in addition to the 7 constant in Eq. flT^ . and 
which is given by the complex argument. The expression of this additional phase is given 
by Eq. ([T]) in the introduction, and is held- and temperature-dependent and proportional to 
the curvature factor k. It is worthwhile to notice that (p{x) changes its sign at x = 1.606115. 
Besides, when x is large, i.e. B/T goes to zero, (p{x) ~ 2 — x. Nevertheless x is much larger 
than 2 in this range, hence (p{x) ~ —x. Therefore, the phase factor is given by 


97^4^2 tj 7^2 

e B 


(18) 


In the opposite case, for very large B/T ratio or very small effective mass, (p{x) ~ x/3, and 
the Onsager phase can be approximated by 


97]-4^2 fc 7^2 

3e ^ ^ B 


(19) 


The above asymptotic expressions Eq. flT 8 |) and Eq. indicate a square temperature de¬ 
pendence and an inverse held dependence of (j)p both in the high an low held ranges. Besides, 
within the tight binding model, 5'"(p) = / (Stt/i^), indicating largest ehect for large unit 

cell and ehective mass. These two features are achieved in organic metals. As an example. 


unit cell area in the cone 


range rria = 3 to 3.51 23 


ucting plane as large as 108 A^j^, 


2211 and ehective masses in the 


27| have been reported for the a orbit of k-(ET) 2 [Cu(NCS )]2 (where 


ET stands for the bis-ethelyne-dithio-tetrathia-fulvalene molecule), k can be estimated from 
band structure calculations and crystallographic data of k-(ET) 2 [Cu(NCS) 2 ] [ 2 l|,l 22 |, yielding 
e.g. K ~ 5 X 10^®m“^J“^ in the M — T direction, relevant to the a orbit. 

As displayed in Fig. [U we have plotted the Onsager phase 0i versus held for several 
temperature values. The phase is very small at high held. Oppositely, large ehects can be 
observed at low held and high temperature. The dashed line determines the threshold limit 
when the Landau level gap is equal to the thermal huctuations, hwe — ksT, and above 
which the oscillations can be observed at higher helds. Additionally, there is another limit 
based on the Dingle temperature To, which also imposes a minimum threshold for the held 
above which the oscillations can be seen. If we take To = 0.5K (the best temperature 
for k-CuNCS [^) one obtains a threshold held equal to 7.5 T when uJcT = 1, and where 
r = h/ksTjj is the scattering time. One however observes oscillations below this threshold 
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FIG. 1. Field dependence of the Onsager phase (l)i at various temperatures deduced from Eq. ([T]) for 
parameters relevant to the a orbit of the quasi-two dimensional organic metals k-(ET) 2 [Cu(NCS) 2 ] 
{d‘^S/dE‘^ = 5 X and m* = 3.2 rUe). The dashed line marks the field values such as 

ksT ~ hwc) which can be regarded as an estimation of the lower limit for the field above which 
oscillations can be observed. 


limit in some compounds, when for example uJcT is equal to 0.28 29|. To our knowledge, 
quantum oscillations in organic metals have been up to now observed above several teslas at 
liquid helium temperatures. Therefore, the studied effect could only be observed provided 
clean compounds with very small scattering rate are synthesized. In this low held range, 
where 0i/27r varies by several units, can be rewritten as = —27i'p5F/B where 6F is 
given by: 


6F 




(20) 
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As a typical example, the tight binding model yields 

6F = {kBTf{Ti‘^/2A)m*^Ao/{en^). 

For m* = rUe and Aq = lA^, i.e. for parameters close to those of elemental metals, 6F = 
3.462 X which remains negligibly small, even at high temperature. For the above 

considered organic metal, temperature as high as few tens of a Kelvin is nonetheless necessary 
to get a frequency shift of few T. Namely, 6F value as low as 1.5 T should be observed at 
20 K which could hardly be detected owing to the oscillation frequency value Fa = 600 T 
yielding 6F/F = 0.25 %. 


III. DIRAC FERMIONS 


In this section, we consider the case of Dirac fermions such as observer 
graphene which has been intensively studied (for a review, see e.g. Refs. 


in monolayer 


30 


3l|). The 


data are in agreement with a linear band dispersion, the curvature parameter being given 


by S''{E) = 27r/(huir)^ with a Fermi velocity vp = 10®ms“^ 3, 32,[33|. Shubnikov-de Haas 
(SdH) oscillations with effective mass in the range 7 x 10“^ rUe to 0.04 rrie, depending on the 
carrier concentration driven by bias voltage, have been observed at temperatures either up 


to 50 K 


or above 100 K 


32|, allowing to expect detection of larger values of the phase 


Contrary to the parabolic case, the Landau level energies En are not linear with index 
n but are given by = Ei^/n^ with Ef = 2ehvpB. Since E'^ oc E'^ oc n~^ and 

E” oc the two terms in the right hand side of Eq. ([5]) are diverging when n goes to 

zero. Therefore, the intermediate step to compute Dose using two integrations by parts in 
Eq. ([5]) cannot be performed here, since these contributions makes the integral in Eq. (l5ll 
divergent. Taking into account the band degeneracy qd = Am. the case of graphene, we can 
write instead 


A 


-2gDbReY^ 

p=i ■ 


exp(2i7rpn) E'^V- 
1 + exp[/3{En — /i)] 2i7ip 



exp{2i'KpE‘^/El) ^ 
1 + exp[f3{E - /i)] 


(21) 
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An integration by parts is possible if we integrate exp(2z7rpii^^/i?^) using the error function 
in the complex plane 




A 


2gDb Re ^ 




TT /■°°er 


i{y/—2i7rpE / El 


^ 4(-2i7rp)3/2 2 Jq cosh^[/3(E -/i)/2] 


E. 


( 22 ) 


We can set as before x = (3{E — p)/2n and obtain 




A 




El 


^3/2 roo erf ^v'-2f7rp{/i + 27ix/I3}/Ei'^ 


X. (23) 


p=i 2(-2i7rp)3/2 2 y_/3^/(2,r) cosh^(7ra:) 

The dominant part of the magnetization is obtained by differentiation of the erf function 
with respect to minus B in the previous expression, since the dominant oscillating term 
comes from this function 


^osc — QD' 




h 


Zp, 


P>1 


Zn — 


TT 


^H-oo 


’ —‘Icrji 


exp{2i7ipE/ B) 
2i7Tp 


Tx\ exp{2i7ip[TX + r^a;^/(4cr)]) 


^ ^ 2a) 


cosh^(7rx) 


X. 


(24) 


Zp is normalized such that Zp = 1 at zero temperature, i.e. r = 0, using the integral 
value TT cosh(7rx)“^x = 1. In this expression, the dominant frequency E is dehned by 
fji^/Ef = E/B = a, and the damping factor r = A7ip/{(]Ef) which is temperature- and 
held-dependent. These parameters are expressed as 


E = 




2ehvp ’ 


r = 


27r/i 


{3ehvpB 


a = 


B'‘ 


T 


a 


dvr 

/3/i' 


(25) 


In reference 3^, the authors obtained the expression for the zero-temperature thermody¬ 
namic potential, given here by Eq. fl23p . in terms of a series involving Bernoulli polynomials. 


up to the fourth power in magnetic field. This is very similar to the expansion obtained for 
the general case from Eqs. ([6]) and flTTl) in the sense that curvature induces an expansion 
over the cyclotronic frequency. As in Eq. (fTTD . Bernoulli polynomials Bn{x) of order n are 
periodic functions of x with harmonics p that decay like p~'^. In Eq. ffTTl) this corresponds 
to an expansion up to n = 3 that includes the temperature dependence. In the case of 
Dirac fermions, this expansion is hidden in the factor Zp whose dependence on the harmon¬ 
ics index p can be explicitly determined by expanding the argument in the exponential 
of Eq. fl2T)) as a series. The remaining integrals can in principle be performed in the com¬ 
plex plane yielding a power series in terms of 1/p. In the limit where a ^ r, we have 
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Zp = npr/ sinh( 7 rpr) = i?p( 7 rr). We may then identify vrr with 2'k ‘^/as in Eq. 
by dehning the cyclotronic frequency Uc = eB/m* and effective relativistic mass m* = pi/v\. 
Otherwise, the amplitude factor Zp = 74 pexp(z 0 p) is complex and possesses a modulus Ap 
and non-zero argument <pp contributing to the global Onsager phase of the oscillations. 

We may apply the residue theorem to the previous integral, in the limit where the lower 
bound —2cT/r is large, or typically T < m*v‘j^/kB- The singularities of the function cosh(a:)“^ 
are still located a.t Xn = i{n + \) and we obtain 


Z ~ 

ZJn — 


^ < 2p7rr 


n>0 


1 + 


iryn -|- 


1 ^ r 1 

1 r 

2ia 

> exp 


— 2'KpT(n + 5 ) — 2i7rp—{n -|- .(26) 



FIG. 2. Field dependence of the Onsager phase (pi for Dirac fermions derived from Eq. ([2]) with 
Fermi velocity vp = 10®ms“^ and effective mass m* = 0.02 at T = 10 K, 50 K and 100 K. 

This expression can be compared to Eq. (1161) provided substitutions k = 1/p and uj = 
/Aa are made. Except for the argument of the exponential terms, Eq. fl26l) is different from 
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the general case, yielding therefore a modihed form of the Onsager phase (j)p and function 
(p{x) dehned in Eq. ffT7|) . given by Eq. ([2]). In this case, the series expansion valid at low 
m*T/B ratio gives (p{x) ~ x instead of x/3, and the phase (pp is inversely proportional to 
the magnetic held. The other limit as m*T/B is large, is as previously negative, although 
(p{x) ~ 4 — a; instead of2 — x. As in the previous case, x is much larger than 4 in this range, 
hence (p{x) ~ —x. Therefore, the phase factor is given by 

2ir^p f ksT^ ^ 


heB 


(¥)■ 


(27) 


where the sign plus and minus corresponds to the high and low held limit, respectively. 
According to Eq. fl27|) . pp depends on the Fermi velocity, instead of the ehective mass and 
band curvature, at both high and low held. Field dependence of (pi is displayed in Fig. [2] 


for experimental parameters relevant to graphene 18 


As reported in the preceding 


section, negligibly small values are obtained in the high held range. In the low held range, 
the same behaviour as for organic metals is obtained. Namely, a signihcant drop down of pi 
at low held and high temperature. Nevertheless, the frequency shift, which can be written 


6F/F = 271^ 


kuT 


, remains small (about 8 % at 100 K in the example of Fig. |2]) although 


higher than in the case of the organic metal considered in the preceding section. 


IV. SUMMARY AND CONCLUSION 


We have evidenced that a non-parabolicity of the band dispersion versus momentum yield 
a held-dependent shift of the Onsager phase of quantum oscillations, i.e. a frequency shift. 
The analytical expressions of this phase shift depend on the nature of the dispersion band, 
either mainly parabolic with corrective terms to the parabolicity (organic conductors) or 
linear (Dirac fermions). It is demonstrated that, in both cases, the effect is amplihed at low 
magnetic held, high temperature and either large effective mass and unit cell area (organic 
conductors) or small Fermi velocity (Dirac fermions). However, the observed effect is small 
within the experimental conditions where quantum oscillations are observed. Largest effect 
are nevertheless observed for Dirac fermions. Indeed, the deviation from parabolicity of 
their band dispersion is the most signihcant. Since the phase shift is expected to be the 


largest for small Fermi ve’ 


in a-(BEDT-TTF) 2 l 3 


35 


ocity, Dirac fermions in organic conductors, such as observed 


36| . with large unit cell, hence small Fermi velocity appear as 
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promising candidates. It is expected that samples should be clean enough to observe the large 
phase shift deviations at low helds. This typically corresponds to the integer quantum Hall 
regime, which requires a proper treatment of the disorder effects but are not that relevant 
for the magnetization oscillations, in contrast to transport coefficients. The temperature- 
dependent phase shifts are not expected to be fundamentally affected. Indeed, the principal 
modification for the magnetization oscillations in the quantum Hall regime will take place 
at the level of the disorder Dingle factor (associated with a Lorentzian distribution of the 
Landau bands) and may be no more valid in low-dimensions, and can be typically replaced by 
a Gaussian reduction factor (with a Gaussian dependence on the harmonics p), which reflects 
a Gaussian shape of the Landau bands as often reported experimentally for the density of 
states in the integral quantum Hall regime [^. Also, it may be important to consider 
theoretical Fermi surfaces with singular behavior which is affected by large curvature and 
therefore may display important phase shifts. 
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